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Abstract 

We present approximate formulas for the tensor BB, EE, TT, and TE multipole coefficients 
for large multipole order I. The error in using the approximate formula for the BB multipole 
coefficients is less than cosmic variance for i > 10. These approximate formulas make various 
qualitative properties of the calculated multipole coefficients transparent: specifically, they 
show that, whatever values are chosen for cosmological parameters, the tensor EE multipole 
coefficients will always be larger than the BB coefficients for all i > 15, and that these 
coefficients will approach each other for i <C 100. These approximations also make clear how 
these multipole coefficients depend on cosmological parameters. 
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I. Introduction 



Tensor fluctuations are a prime target for future observations of the cosmic microwave 
background, because if detected they can provide a conclusive verification of the theory of 
inflation and a unique tool for exploring the details of this theory. The contribution of these 
fluctuations to the correlations of temperature and polarization correlations is well known. 
They have the multipole coefficients^ 
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Our notation here is consistent with that of reference j2j. To be explicit: To is the microwave 
background temperature at the present time t ; P(t) = uj c (t) exp[— j/° u c (t') dt'} is the prob- 
ability distribution of last scattering, with u c (t) the photon collision frequency; t\ is any 
time taken early enough before recombination so that any photon present at t\ would have 
collided many times before the present; r[t) = J* t *° dt'/a(t') is the co-moving radial coordinate 



1 These formulas are equivalent to those of Zaldarriaga and Seljak pQ. Their gravitational wave amplitude 
h and power spectral function Ph{k) are related to our gravitational wave amplitude D q (t) by h\fP~h = T>/2. 
In consequence, their function ty\fP~h is 1/4 times our source function ty. 



of a source from which light emitted at time t would reach us at the origin at the present time 
to; and ty(q,t) is the "source function," which is customarily calculated from a hierarchy of 
equations for partial-wave amplitudes: 

A$(g,f) + a(2 / +1) ((* + - ^ { T,L(Q,t)) 

= (-2P 9 (t) +w c (t)tt(g,*)) 5 lfi -uj c (t)A { p E (q,t) , (5) 

A£J(g,t) + q(2 / + i) ((' + *) A wfa ~ *)) 

= -LU c (t)V(q,t) 5 efi ~ uJ c (t)A ( p}(q,t) , (6) 

with 

*(?,*) = 7^ A SM) + ^A$(g,t) + ^A£](g,t) - -A^g,*) 

+ ^Ag(g,t)- ^Ag(g,t). (7) 

Here T> q (t) is the gravitational wave amplitude (apart from terms that decay outside the 
horizon), defined by 

<fcfc(x,t) = £ / d 3 ge 4q - x /3(q,±2)e y (g,±2)D g (t) , (8) 

with /3(q, ±2) and ey(g, ±2) the stochastic parameter and polarization tensor for helicity 
±2, normalized so that 

(/3(q,A)^(q',A'))=5 AA ^ 3 (q-q'), (9) 

and for q in the 3-direction 

en(g,±2) = -e 22 (g,±2) = 1/^2, e 12 (g, ±2) = e 21 (g, ±2) = ±i/V2 . (10) 
Finally, d(g, t) is the quantity 



d(g, t) = exp 



jT t0 rft'^(t')] (i>,(t) - ^ c (o*(g,*)) . (n; 



Aside from the treatment of the tensor mode as a first-order perturbation, and the as- 
sumption of purely elastic Thomson scattering, Eqs. (pQ)-(jl]) may be regarded as exact. 



They serve as the basis of computer programs such as CMBfast and CAMB, that are used 
to compare observations of microwave background polarization and temperature fluctuations 
with models that predict values for the gravitational wave amplitude T> q (t). But they are 
not very transparent. 

For one thing, as shown in Figure [fl computer calculations using Eqs. §B) and (T5]) yield 
results for C EE e and C BB e that are of the same order of magnitude, and nearly equal for 
£ < 100, while C T EEl > C T BBt for all t > 15. 
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Figure 1: Comparison of C EE e and C BB e , in (/uK) 2 . 

Of course computer calculations can only show this for specific choices of cosmological 
parameters. (The cosmological parameters used in Figure 1 are described below.) It would 
be impossible to conclude just by inspection of Eqs. (CD) and (j2D that these are general prop- 
erties of the multipole coefficients, independent of the choice of cosmological parameters. In 
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this paper we present successive approximations that make these properties apparent, and 
that, at the cost of only a small additional loss in accuracy, also clarify how the multipole 
coefficients depend on various cosmological parameters. 



II. The Large-£ Approximation 



We can approximate Eqs. (JH) — by much simpler and more transparent formulas, by 
using an asymptotic formula |3J for the spherical Bessel functions o 
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where v = I + 1/2, and cos b = u/p, with < b < n/2. This approximation is valid for 
\v 2 — p 2 \ ^> z/ 4 / 3 . Hence for £ ^> 1, this formula can be used over most of the ranges of 



integration in Eqs. (JTJ) — (HI) • Furthermore, for p > v ^> 1 the phase z/(tan6 — b) in Eq. ( TT2l 
is a very rapidly increasing function of p, so the derivatives in Eqs. ([I])-© can be taken to 
act chiefly on this phase: 
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2 After the preprint of this work was first circulated we learned that the same approximation is used 
by J. R. Pritchard and M. Kamionkowski, Ann. Phys. 318, 2 (2005) [astro-ph/0412581| . However, after 
making this approximation they make further approximations that are quite different from ours, and that 
lead to a divergence in the integral over wave number, which must be dealt with by an arbitrary cut-off. 
The error introduced by their approximation is comparable to the one introduced by our last approximation 
given by Eqs. (|23|) and (|24| . Another approximation that consists in averaging over the rapid oscillations in 
the square of the Bessel functions leading to results similar to our last approximation was proposed by M. 



Zaldarriaga and D. D. Harari, Phys. Rev. D 52, 3276 (1995) astro-ph/9504085|. An analytic expression 
for the contribution of the tensor modes to the temperature multipole coefficients approximately valid for 
1 I < 50 obtained using a similar average was given by A. A. Starobinsky, Sov. Astron. Lett. 11, 133 
(1985) 
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Then Eqs. ©-® become, for z/ = £ + 1/2 > 1, 
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In evaluating both the exact and approximate expressions, instead of calculating the 
source function ty(q,t) by truncating the Boltzmann hierarchy (JSJ), ([6]), we use the integral 
equation [2], H] 
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where K(y) and F(v) are the functions 

K(v) = j 2 (v)/v 2 , F(v) = j (v) - 2 3l (v)/v + 2j 2 (v)/v 2 . (20) 

(This is not an approximation; in principle it should give the same results as the truncated 
Boltzmann hierarchy used by CMBfast and CAMB, aside from the supposedly small errors 
produced by the truncation. In fact, our method gives results that differ by a few percent 
from both CMBfast and CAMB, but CMBfast and CAMB give results for both C^ El and 
Cg B e that differ by similar amounts from each other, especially for large t. At this point we 
are not able to tell which of the three methods is the most reliable.) The specific cosmological 
model chosen for this and all other numerical calculations in this paper is consistent with 
current observations: We assume zero spatial curvature and constant vacuum energy, with 
density parameters for baryons, cold dark matter, and dark energy given by Qb = 0.0432, 
Qcdm — 0.213, £7a = 0.743. We take the reduced Hubble constant as h = 0.72 and 
the present microwave background temperature as T = 2.725K. In calculating the photon 
collsion frequency, we use the recfast recombination code [8], with helium abundance Y = 
0.24. The gravitational field amplitude outside the horizon is taken as 

\V q \ 2 = 4.34 x 10- n <r 3 

corresponding to tit = 0, A = 0.739, and a tensor/scalar ratio r = 1 in the notation of [9|. 
(The values of the parameters used are the rounded maximum likelihood values from the 
full N-dimensional likelihood analysis and may differ slightly from the marginalized values 
quoted in [9]). The gravitational wave amplitude T> q (t) is calculated including the damping 
due to neutrino anisotropic inertia, as in [10J; the effects of photon anisotropic inertia are 
negligible. Reionization is ignored. To take into account a finite optical depth r of the 
reionized plasma or a different value of r, for I > 10 it is only necessary to multiply the 
multipole coefficients given here by rexp(— 2r). The approximate results obtained in this 
way from Eqs. (|15|) - (|18p are compared with the exact formulas (JTJ-(jlj) in Figures [2HSJ 

We can gain further simplicity and transparency in the formulas for the EE and BB mul- 
tipole coefficients by using another approximation that actually leads to improved accuracy 
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for C BB e . The last-scattering probability distribution P(t) is concentrated around a time 
tL, corresponding to a redshift zl — 1090. For any q of the same order of magnitude as 
zVK^i)) the quantity b = cos -1 {y/qr{t)^ does not vary appreciably for t within the range 
in which P(t) is appreciable. Hence we can set r(t) equal to everywhere except in the 
phase u(ta.nb — b), which for v ^> 1 does vary over a wide range in this interval. Fur- 
thermore, because z/(tan6 — b) varies over a wide range for v 3> 1, the difference between 
cos[z/(tan6 — b) — 7r/4] and sin[i/(tan6 — b) — 7r/4] is immaterial, and we can replace both 
with cos[z/(tan& — 6)]. Making these replacements in Eqs. (fT5|) and (fl6|) gives 
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(We have not set 6 = b^ in the factors l/vsin6 in both integrals over t, in order to 
avoid a divergence in the integration over q at q = vjr^. This factor does not introduce a 
divergence in the integrals over time, because dt oc sinbdb.) 

These approximate formulas are compared with results of the exact formulas ([1]) and ([2]) 
in Figures [6] and [7J The approximate result (|22|) for C BB e agrees with the exact result (T5]) 
to about 1% for all £ > 10, which is better than cosmic variance. The approximate result 
(|2~TT) for C BE e is not quite as accurate; it agrees with the exact result (CQ) to better than 
about 14% for all £ > 10. These approximations are evidently accurate enough for us to 
draw qualitative conclusions about the EE and BB multipole coefficients. 

One immediate consequence is that, since (1 + sin 2 6^) 2 > 4sin 2 6^ for all real b^, we 
expect that C EE e > C BB t for all £ large enough to justify our approximations. Also, since 
\&(g, ti) falls off for wave lengths that come into the horizon before matter-radiation equality, 
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we expect that for relatively small i (say, I < 100) the integrals over q are dominated by 
values for which cos 6/, is small, so that (1 + sin 2 ^) 2 ~ 4sin 2 &r J , and hence C% Ei — Cg B e 
for such £. As mentioned in Section I, and shown in Figure [H, both properties are observed 
in the output of numerical calculations based on the accurate formulas ([!]) and (J2J). 



III. Parameter-Dependence of the EE and BB Correlations 

With one further approximation, we can find reasonably accurate formulas for C EE £ and 
Cg B e that reveal the way that these coefficients depend on various cosmological parameters. 
We write the squared time integrals in Eqs. ff2T]) and (1221) as double integrals over times t 
and t' and write 
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where cosfr = vjqrit) and cos6' = vjqr(t'\ For v ^> 1, and qr(t) and qr(t') both of 
order z/, the second term on the right oscillates very rapidly, and hence may be neglected in 
the integral over t and t'. On the other hand, because P(t) and P{t') are sharply peaked 
around the same time ti, the argument of the first cosine on the right is small where P(t) 
and P(t') are appreciable, so this cosine may be replaced with unity. Then (now dropping 
the distinction between v and £), Eqs. ( 1211 and ( 1221) become 
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This approximation is compared with the results of the exact formulas ([I]) and (|2J) in 
Figures [8] and [91 As shown there, the fractional error here is less than about 20% for 
10 < £ < 600, but it becomes larger for larger values of £, where the multipole coefficients 
become quite small. 

Eqs. fl23l) and (j2"4"|) are useful in revealing the parameter dependence of these multipole 
coefficients. Where the last-scattering probability distribution P(t) is appreciable, the only 
cosmological parameters on which either P(t) or the source function ^(q, t) depend are the 
baryonic and matter density parameters fish 2 and Q]\,fh 2 , as well as the present microwave 
background temperature To. All dependence of the multipole coefficients on Hq or the 
curvature flxh 2 or the vacuum energy Q^h 2 is contained in the function r(t). But Eqs. (1231) 
and (124"|) show that r(t) and £ enter in the multipole coefficients only in the combination 
r(t)/£. Hence, with f^/i 2 , ^Aih 2 , and T fixed, to a good approximation C EE e and C EBl 
depend on H , Qxh 2 , and f^^- 2 only through their effect on the scale of the ^-dependence 
of C EE e and C BB e . Furthermore, since P(t) is sharply peaked at the time of last scattering, 
just as for scalar modes there is a high degree of degeneracy here: for £ > 10 the coefficients 
C EEi and C BBl depend on H , Q K h 2 , and ^a^ 2 only through a single parameter, the radius 
r{ti) of the surface of last acatttering. Of course, the degeneracy here is not as important 
as it is for scalar modes, because tensor modes when discovered will be studied primarily 
for the purpose of measuring the tensor/scalar ratio r and the tensor slope tit, rather than 
other cosmological parameters. 

We are grateful to Eiichiro Komatsu for frequent helpful conversations. This material 
is based upon work supported by the National Science Foundation under Grant No. PHY- 
0455649 and with support from The Robert A. Welch Foundation, Grant No. F-0014. 
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Figure 2: Comparison of formulas for C EE E . The solid line is the result of using the exact 
expression (pQ); the dashed line is the result of using the approximation ffl5l) . Figures 2-5 show 
the degree of accuracy of the large-£ approximation by itself, without further approximations. 
In this and all other figures, all calculations are done using the cosmological parameters given 
in Section II, and the units of the vertical axis are square microKelvins. 
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Figure 5: Comparison of formulas for C^ Te . The solid line is the result of using the exact 
expression (TjJ; the dashed line is the result of using the approximation ffT8l) . 
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Figure 6: Comparison of formulas for C EE t . The solid line is the result of using the exact 
expression (CD); the dashed line is the result of using the approximation ( l2~li) . Figures 6 and 7 
show the degree of accuracy of the combined approximations that we use to show analytically 
that C EEe > C] 3B l . 
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Figure 7: Comparison of formulas for C BB t . The solid line is the result of using the exact 
expression (J2J); the dashed line is the result of using the approximation (1221) . This is our best 
approximation for C BBe 
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Figure 8: Comparison of formulas for C EE t . The solid line is the result of using the exact 
expression (CQ); the dashed line is the result of using the approximation ( |23|) . This figure 
shows the degree of accuracy of the further approximations used to explore the parameter- 
dependence of C EE t 
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Figure 9: Comparison of formulas for Cj 3B e . The solid line is the result of using the exact 
expression (j2J); the dashed line is the result of using the approximation (1231) . This figure 
shows the degree of accuracy of the further approximations used to explore the parameter- 
dependence of Cg B t 
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